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Abstract 



We present numerical solutions to the extended Doering-Constantin vari- 
ational principle for upper bounds on the energy dissipation rate in turbu- 
lent plane Couette flow. Using the compound matrix technique in order to 
reformulate this principle's spectral constraint, we derive a system of equa- 
tions that is amenable to numerical treatment in the entire range from low 
to asymptotically high Reynolds numbers. Our variational bound exhibits a 
minimum at intermediate Reynolds numbers, and reproduces the Busse bound 
in the asymptotic regime. As a consequence of a bifurcation of the minimizing 
wavenumbers, there exist two length scales that determine the optimal upper 
bound: the effective width of the variational profile's boundary segments, and 
the extension of their flat interior part. 
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I. INTRODUCTION 



Can one derive upper bounds on the rate of energy dissipation in turbulent shear flows 
that are both mathematically rigorous and physically meaningful, i.e. "sharp"? This ques- 
tion is of enormous conceptual importance in the theory of turbulence, since the energy 
dissipation rate is the central object in classical turbulence theories and its behavior 

at high Reynolds numbers may reveal whether there are intermittency corrections to clas- 
sical scaling Since full numerical simulations of shear flows with Reynolds numbers 
of the order of 10 6 or above are out of reach, and will remain so in the foreseeable future, 
rigorous estimates of the dissipation rate derived directly from the Navier-Stokes equations 
are one of the few tools left to the theorist for approaching such questions. 

The idea of bounding the rate of energy dissipation with the help of rigorous inequal- 
ities, rather than trying to solve the equations of motion under plausible, but essentially 
uncontrolled approximations, saw its first prime in the early seventies, with the formula- 
tion and further exploration of the Howard-Busse theory || . In particular, Busse [§-[3 \ 
was able to derive an asymptotic upper bound on energy dissipation in plane Couette flow 
that is in accordance with classical theories, but lies roughly an order of magnitude above 
experimentally measured data [fl~2|Jl~3[1 . 

After this state of affairs had remained unchanged for about 25 years, Doering and 
Constantin put forward a new variational principle for computing quantities characterizing 



turbulent flows flI^|T6|j. Whereas the first applications of this principle to the plane Couette 



problem [14,17|)lq1 did not exhaust the principle and produced bounds on the dissipation 



rate that could not compete with Busse's, subsequent analysis by Kerswell |T^] revealed 
that the Doering-Constantin principle, after being extended along the lines worked out in 
Ref. pi, actually reproduces the Busse bound in the asymptotic regime. Kerswell arrived 



at this conclusion by mapping the extended Doering-Constantin principle to another vari- 
ational principle that also occurs in Busse's theory, and then applying Busse's method of 



solution. However, as emphasized by Kerswell himself |19[ , there is still no mathematical 
proof that Busse's solution is really extremal, although this is generally believed to be the 
case. Moreover, the Busse bound pertains to the limit of asymptotically high Reynolds num- 
bers only. Hence, what is needed is a solution to the extended Doering-Constantin principle 
that (i) fully exhausts this principle, (ii) spans the entire range from low to asymptotically 
high Reynolds numbers, and (iii) is logically independent from Busse's solution. 

In this paper, we present a full-fledged numerical analysis of the improved Doering- 
Constantin variational principle for the plane Couette flow that meets these requirements. 
The goal here is not merely to derive the best possible bound, but we also wish to explore the 
working principles of the new method. Even if statements concerning the existence of non- 
classical scaling remain out of reach at present, the potential power of the new variational 
principle lies in its great conceptual clarity. Once thoroughly understood, it might serve as 
a germ for the formulation of a still more advanced principle that could actually allow one 
to attack the scaling questions. 

We organize our material as follows: after having stated the variational principle in 
Sec. [TIJ Sec. |TIT] will be devoted to the resolution of the most demanding technical difficulty, 
the implementation of the principle's spectral constraint. It is this constraint that seems 
to have hindered previous numerical work [^TJ to reach the asymptotic regime. Our results 
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will then be described in Sec. [TV]; the final Sec. [V| offers a concluding discussion. Some 
important technical details concerning the derivation of the system of differential equations 
actually used in the numerical computations can be found in Appendix [A|, while Appendix [FJ 
contains a brief, but instructive example that illustrates the key points of our method, and 
underlines its accuracy, without too much technical burden. 

Since the numerical work is fairly involved, it would hardly have been possible to carry 
through the entire investigation without some analytical guidance. Our guide is the example 
problem of the Couette flow without spanwise degrees of freedom, referred to as the restricted 
Couette problem in the following, which has been treated in detail in Ref. |[22|| , and which 
shows several features that will reappear here in the context of the unrestricted problem. 



II. THE VARIATIONAL PRINCIPLE 



We consider a layer of an incompressible fluid with kinematic viscosity v that is confined 
between two infinitely extended, parallel rigid plates separated by a distance h. The lower 
plate, coinciding with the plane z = of a Cartesian coordinate system, is fixed, whereas 
the upper plate at z — h is sheared with constant velocity U in the positive x-direction. The 
dynamics of the fluid's velocity field u(x, t) are determined by the equations 



d t u + u ■ Vu + Vp = v Au (Navier-Stokes equations) , 
V • u = (incompressibility) , 
u(x, y, 0, t) = 0, u(x, y, h,t) = £/x (no-slip boundary conditions); 



(1) 
(2) 
(3) 



here p denotes the kinematic pressure and x is the unit vector in x-direction. In addition, in 
both x- and y-direction periodic boundary conditions (b.c.) are imposed on u and p. The 
time-averaged rate of dissipated energy per mass is given by 



e T 



7 



dt 



v 
ft 



d x 



h3=x,V,z 



(4) 



where Q is the periodicity volume. Our aim is to derive a mathematically rigorous upper 
bound on the long-time limit e of et, 



e = lim Et. 

Equivalently, we consider the non-dimensionalized dissipation rate 

£ 



c £ (Re) 



U^h- 1 



(5) 



(6) 



where Re = Uh/u is the Reynolds number. For calculating the bound we will employ the 
background flow method as recently put forward by Doering and Constantin [|14} -|I7|, and 
developed further in Ref. |20| . 

It is known [14[] that the rate of energy dissipated by the laminar flow u(x, t) = (Uz/K)5t 
provides a rigorous lower bound on c e for all Re, namely c e {Re) > £e(Re) = Re -1 . 
Moreover, energy stability theory |23|j2^1 guarantees that the laminar flow is the only possible 
stable flow at least up to the energy stability limit Rces- Hence, we have the exact identity 
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c £ (Re 



Re 



for Re < Re ES ~ 82.65. 



(7) 



From an early generalization of energy stability theory, the Optimum Theory devised by 
Busse pHTlj, one obtains an approximate asymptotic upper bound [p5|1 : 



lim 

Re^o 



c £ {Re) < 0.010. 



(8) 



Note that the right hand side of this inequality does not depend on the Reynolds number, so 
that this bound is in accordance with classical turbulence theories [0^0] • (For a discussion of 
the connection between the .Re-dependence of c e {Re ) and possible intermittency corrections, 
see Ref. [§.) 

The background flow approach can likewise be considered as a generalization of energy 
stability theory, but it has a distinctly different twist. Instead of decomposing the velocity 
field u into a ^-dependent, plane-averaged mean flow and the fluctuations around this flow, 
as Busse did in his Optimum Theory, Doering and Constantin [T7J take up an idea by 
Hopf p6| and write 



u(x,t) = U(x) + v(x,t) 



(9) 



Now the stationary and divergence- free auxiliary field U(x), dubbed the background flow, 
has to carry the physical b.c, namely U(x, y, 0) = 0, U(x, y, h) = £/x, and U(x) is periodic 
in x- and y-direction. But apart from these fairly mild specifications, the background flow 
remains completely arbitrary. Hence, this decomposition (^) serves as the starting point 
of a variational principle [14|: plug it into the expression (||) for the energy dissipation 
rate, get rid of the deviations v(x, t) with the help of both the equations of motion (|l])-([3]) 
and rigorous inequalities that bound the dissipation rate from above, and then adjust the 
background flow such that this upper bound becomes as low as possible. 

For carrying through this program, we restrict ourselves to background flows that are 
given by a merely height-dependent profile 0(C), which should reflect the symmetry of the 
Couette geometry. That is, we admit only background flows of the form 



U(x) = £ty(C)i, 



(10) 



with 



0(0) = 0, 0(1) = i, 0(O = i-0(i-C), (ii) 

where ( = z/h is the dimensionless coordinate in cross-stream direction. The resulting 



inequality that bounds the dissipation rate then reads |20 

,2 



c £ {Re) < 



1 + I7^T) fl W 



Re~\ 



(12) 



where a > 1 is a dimensionless balance parameter, and D{(p} denotes the profile functional 



± 

d{0} = fdcmo? 



i. 



(13) 
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The bound on the dissipation rate provided by the r.h.s. of the inequality (|12D can be 
minimized by varying both the profile <fi and the balance parameter a. However, there is 
a substantial technical difficulty: the inequality ( |7I| j is valid only as long as the profile <fi 
satisfies a spectral constraint Hl4| . This spectral constraint demands that all eigenvalues A 



of the linear eigenvalue problem 

/0 1 

AV = -2 h 2 AV + R<j>' | V + VP, 

V 1 

o = v-v, v(on) = o (14) 

for the stationary velocity fields V be positive. Since the background flow already carries the 
physical boundary conditions, V has to satisfy homogeneous b.c, that is, V has to vanish 
on the material walls dQ at z = and z = h, and to be periodic in x- and ^-direction. 
Note that the balance parameter a enters into this eigenvalue problem through the rescaled 
Reynolds number 

R = -^—Re. (15) 



In effect, the spectral constraint (|TJ) restricts, for each Reynolds number, the combina- 
tions of profiles and balance parameters that are admissible to the variational principle (|T^|) . 
The highly non-trivial task of evaluating this constraint consists in determining, for each can- 
didate profile separately, that rescaled Reynolds number R where the smallest eigenvalue 
passes through zero. We will denote this number as R c {4>} in the following. 

Provided we know how to compute R c {4>} — this task will be taken up in the following 
section — , each leads to a rigorous upper bound on c e in the interval < Re < R c {<p}. 
After optimization of the balance parameter, this bounds adopts the form |20| 



C £ < 



D{(j)}} Re" 1 for < Re < ±P C {0} 

(16) 



1 , D{<t>}R c {<t>} 



4(R c {4>}-Re)Re 



Re' 1 for |P C {0} < Re < R c {<j)} 



Figure [I] shows the graph of such an upper bound produced by a generic profile <fi. When 
a suitable class of test profiles has been selected, the lower envelope of all graphs stemming 
from the individual profiles then is the optimal upper bound on c E (Re ) that can be obtained 
from this class, i.e., the solution to the variational principle ( |T2"D within this class. 



III. EVALUATION OF THE SPECTRAL CONSTRAINT 

We non-dimensionalize the problem by choosing the gapwidth h as the unit of length. 
By virtue of the linearity of the equations (|i~4]), and the homogeneous b.c. imposed on the 
eigenvectors V, it is not necessary to specify a scale of velocity explicitly. For the sake of 
notational simplicity we denote all non-dimensionalized quantities by the same symbols as 
their dimension-carrying counterparts. For instance, the velocity field V(x) now is regarded 
as a dimensionless vector function of the dimensionless coordinates x, y and z. 

Utilizing the periodic b.c, we start from the ansatz 
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V(x) = v(z) e ^ x+k y y \ P(x) ee p(z) e i ^ x+k ^ 
and transform the eigenvalue equations (113) into the system 

Xv x = -2 (d 2 z - k 2 ) v x + R<p'v z + i k x p, 

\v 7 . 



-2 (d 2 - k 2 j v y + i k y p, 
-2 (d 2 z -k 2 )v z + R^'v x + P ', 



= \ k x v x + ik y v y + v' z , 



(17) 
(18) 
(19) 
(20) 



with 



k EE Jkl + k 2 y . 



The corresponding b.c. for the function v(z) read 

v(0) = v(l) = 0. 



(21) 



If k y = we encounter the model problem of a Couette flow without spanwise degrees of 



freedom. Since this restricted problem has already been studied in detail in Ref. [£2 
only need to consider k y ^ here. 
Defining the linear operator 

C ee 2 {d 2 - k 2 ) + A, 
Eqs. flT7|) and flT8|) together with Eq. (JZOj) yield 

p — [£ v 'z + ik x R<p'v z } . 



wc 



Hence, Eq. (|17D can be brought into the form 

1 



k 2 



k 2 R <f>'v z + i k x Cv' z 



(22) 



(23) 



and we finally obtain from Eq. (|T9|) a rather intricate equation for v z {z): 



C 2 (d 2 - k' 



C (d 2 z - k 2 ) v z 



, , (2) 



+ (fc yJ R0') ^ + i^^<2£^ + ^— £^ + 4 



0(3) /0(a)- 



+ 2 



0(4) ^ 0(3)0(2) 



-3 



+ 2 



0, (24) 



where we rely on the assumption that <f)'(z) does not vanish for all z in the interval [0, 1]. 
The corresponding b.c. for v z are given by 
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v z (z) U= ,i = 0, U =0 ,! = 0, and £ (d 2 z - k 2 ) v z (z) \ Z=0A = 0. (25) 



Equations (|2~4 ) and (|25|) constitute a sixth-order Sturm-Liouville eigenvalue problem. If one 
succeeds in solving this problem, that is, in computing an eigenvalue A and its eigenfunction 
v^'(z) for some fixed profile 0, wavenumbers k x and k y , and a rescaled Reynolds number 
R, then one can determine the associated components v^'(z) and v^ x \z), together with 
p( x >(z), from Eqs. (p3|), fll8|), and (p2|), so that an eigenvalue and a corresponding eigenvector 
(x) , P^(x)J of the original eigenvalue problem ( |14|) have been found. 
Since we cannot hope to find an analytical solution to this formidable problem for each 
candidate profile 0, we have to resort to a numerical treatment. This requires to bring the 
Sturm-Liouville problem into a form that can actually be dealt with on a computer. We 
defer the technical details of this reformulation, that rely heavily on the so-called compound 
matrix method [27], to Appendix |A|, and present only the final equations here. Although 



( V «( 3 



these equations may appear, on a cursory glance, even more cumbersome than the original 
problem posed by Eqs. (]24]) and fl25|) , they have actually been found well suited for computing 
the "critical" rescaled Reynolds numbers R c {4>} for all situations of practical interest, as will 
be amply demonstrated in the following section. 
To begin with, we define the abbreviations 



and the auxiliary functions 

(j>M{z) 



<P2{Z) 



2k 

(2) ( 



z 



<P'(z) ' 



h{z) = 3k+ k^izy-k[-Wizy) ' 



p r w it 2 ^(z) ' 



3 



^) s _0U (z) ___^u + 2fc8 ^ ^ ^ 

, x , 1 (3) W 2 /0 (2) (^)\ 2 1 , „ \ \ 2 



2 ' 



(27) 



The compound matrix technique outlined in Appendix |A] then yields the following system 
of 20 ordinary, first-order complex differential equations: 

y[ = k —3 2/1 + 2/2 — 1 2/4 — 1 2/6 
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k 



—3 2/2 + 2/3 + 2/5 — ^ 2/8 



y' 3 = k [-3 2/3 + 2/4 + 2/e] , 



y' 5 



y\ 



2/3-3 2/4 + 2/7 + 1 2/io] + A [-2 2/2 + 2/4] - A2/1 - i Rx 
-3 2/5 + 2/6 + 1 2/9 + 2/ii] , 
■32/6 + 2/7 + 2/8 + 1 2/10 + 2/12 



>l2/2 + 02 (2/1 + I 2/3) 



1 1/6 



3 2/7 + 2/9 + 2/13 + A [-2/1 — 22/5 + 2/7] — i ^ 



^12/5 + o 022/6 + C/12/1 



y's = k [-3 2/8 + 2/9 + 2/u] 



2/9 



A: 



§2/8 



3 2/9 + 2/10 + 2/15] + A [-2/2 + 2/9] + /22/s - i Rx [02 (-2/5 + 1 2/s) + 5-12/2 



y'10 = k [-3 2/10 + 2/ie] + /1 [-2/3 + 2 2/8 + 2/io] + A2/6 - i Rx [-0i2/s - 022/6 + 9m] 



2/u 

2/^2 

2/13 



-3 2/n +2/12 + 1 2/15 + 1 2/17 
-3 2/12 + 2/13 + 2/14 + 1 2/ieJ , 
1 2/12 - 3 2/13 + 2/15 - 1 2/19I + A [-2 2/u + 2/13] + ^_R y 2/i 



+ i R x [-012/H - I 022/12 + #22/1] , 
2/m = k [-3 2/14 + 2/15 + 2/17] , 

2/u - 3 2/15 + 2/ie + 2/18 - \ 2/20J + A2/15 + A2/11 + ^R y 2/2 
(2/11 - \ 2/u) + c/22/2 , 



24 



= fc 
+ 



2/17 

2/18 



2/ie = fc [-3 2/ie + 2/19] + A [2 2/14 + 2/ie] + A2/12 + ^2/3 + i [0i2/i4 + 022/12 + c/22/3] 
& [-3 2/17 + 2/18 - 1 2/20 

fc [1 2/17 - 3 2/18 + 2/19] + A [2/11 + 2/is] + ^2/5 + i Rx [-| 022/17 + (7i2/ii + #22/5 

2/19 = fe [-3 2/19 + 2/2o] + A [2/12 + 2 2/17 + 2/19] + ^R H 2/6 + 1 Rx [012/17 + 512/12 + 522/6] , 



2/20 



-3/C2/20 + A [2/14 + 2/20] - A2/17 + ^2/8 + i Rx [-022/17 + 5i2/i4 + 522/s] 



(28) 



These equations have to be supplemented with the initial conditions required for their nu- 
merical integration from z = to z — 1: 



2/2o(0) = 1, 2/*(0)=0 for 2 = 1,..., 19. 
In addition, we merely have to satisfy the fairly simple boundary condition 

2/i(l) = 0. 



(29) 



(30) 



What matters here is precisely the replacement of the rather troublesome, symmet- 
ric boundary conditions (|25|) accompanying the original Sturm-Liouville problem by the 
asymmetric boundary conditions (^) and (|30|) , since these new boundary conditions lend 
themselves to a standard shooting method: just integrate the system (^) with the initial 
conditions (|29|), monitor the resulting value of 2/1 (1), and adjust the system's parameters 
such that this value becomes zero. 
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The price to pay for this simplification is a rather complicated system of equations. But 
still, the spirit underlying these equations is simple. Since we are not interested in the full 
spectrum of the eigenvalue problem (|T4]) (or, equivalently, in that of the Sturm-Liouville 
problem [25D), but only in that value R c {4>} of R where the lowest eigenvalue associated 
with some profile passes through zero, we were entitled to set A = in the derivation of 
the system (0). We then keep both the profile <j) and the wavevector k fixed when solving 
the initial value problem (p8| , p9|) , and adjust only the parameter R. The smallest value of 
R where yi{l) vanishes then equals that value of R where the smallest eigenvalue A becomes 
zero, for this particular <p and k. We denote this R as 

i?oM(k) • (31) 

The desired number R c {<fr} then is identified as the global minimum, taken over all wavevec- 
tors. Since the wavenumbers k x and k y effectively enter into the system (|28| ) only through 
their squares k 2 x and ky, we have 

i? c {0}= min {i? o {0}(k)}. (32) 

Thus, the evaluation of the spectral constraint, i.e., the determination of the maximal 
rescaled Reynolds number R c {4>} up to which a given profile <ft remains an admissible test 
profile for the variational principle fll2l), requires some labor. One first has to determine, for 
that cf) and a representative set of wavevectors k, the first zero of |/i(l) that occurs when R 
is increased from zero to positive numbers. This yields the values ([31]). In a second step one 
has to locate the global minimum of all these values, which equals R c {4>}. 



IV. RESULTS 



Now that we are in a position to evaluate the spectral constraint (|14]) for an arbitrary 
profile 4>{z), we have to specify a class of test profiles that is likely to exhaust the variational 
principle (0, |16|). Based on the analysis of the restricted Couette problem, we propose the 
following class of variational profiles: 

f Ul-p)+pz-\{l-p){l-z/8) n for 0<z<5 

(j)(z) = I Ul-p)+pz for 5<z<l-5. (33) 

{ \ {^-p)+Pz + \{l-p){l-{l- z)/5) n for \-8<z<\ 

These profiles contain three independent parameters: 

i. ) the boundary layer thickness 8 (0 < 8 < |), 

ii. ) the slope p of the profile in the interior (0 < p < 1), and 

iii. ) the polynomial order n (n = 4, 5, . . .) of the boundary layer segments. 

By construction, each (p is n — 1 times continuously differentiable at the matching points 
z = 5 and z = 1 — 8. Since the right hand side of the system ( PB| ) of differential equations 
contains derivatives of the profile up to the fourth order, the smallest value of n allowed 
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by this system is four. In the case n = 4 one encounters a discontinuity of the fourth 
derivative. But as long as this discontinuity remains finite, it poses no problem for the 
numerical integration. For practical computations we restrict the parameter n to values less 
than or equal to 1000. 

At this point it needs to be emphasized that the variational parameter n plays an im- 
portant role for the numerical solution of the variational principle. If this parameter was 
artificially kept fixed, that is, if one considered the profiles fl33|) with a constant param- 
eter n = 4, say, then one would inevitably be confronted with serious divergences which 
are hard to deal with numerically, and therefore would drastically reduce the maximal at- 
tainable Reynolds numbers. The formal reasons for these difficulties lie in the expressions 
cf)( 2 \z) /(j>'(z) and (f>^(z) /(f>'(z), which enter into the system fl2E| ) through the auxiliary func- 
tions (0). Namely, for high Reynolds numbers the optimized profiles tend to become flat 
in the interior, so that one has to face small denominators. 

Since the solution to this technical problem, achieved with the help of the variational 
parameter n, is of outstanding practical importance, we illustrate it with an example. Let 
us first note that the parameters 5 and n are related to each other by the profile's slope at 
z = (or, equivalently, at z = 1): 

(f)'(0)=p+-(l-p). (34) 

Secondly, we anticipate the high- .Re-scaling behavior of the optimized profiles: we have 

0'(O) ~ aRe and <//(§) = p ~ (3Re~ l 

with constants a and (3, as will be shown later on (cf. Eqs. (|39|) and (|40~D). This implies 
that the ratio n/(25) in Eq. (Q) necessarily has to increase about linearly with Re when 
Re becomes large. We now focus on two different parameter combinations: (a) 5 = 0.5 and 
n = 1000, and (b) 5 = 0.002 and n — 4. In both cases the ratio n/(28) takes on the value 
1000, and what matters now is the way the divergences of ^ 2 \z) /(j)'(z) and (p^(z) /(f)' (z) 
manifest themselves when the profile slope p approaches zero. Case (a) is actually met 
in high- .Re-solutions of the variational principle when the shape of the profiles' boundary 
segments is allowed to vary, whereas this shape is artificially kept fixed in case (b). In Figs. |2| 
and |3| we demonstrate for these two cases the response of (f)^ 2 \z) /(f)' (z) to the variation of 
p from 0.1 to smaller values, depicting only the relevant intervals of z. In case (a) the 
decrease of p merely entails a shift of the smooth ascend from the (negative) value at z = 
to zero. In marked contrast, in case (b) the same decrease of p produces a rather sharp 
singularity in the vicinity of z — 5. The expression — <p^\z) /(f)' (z) shows, qualitatively, the 
same behavior, but even more pronounced. While the functions shown in Fig. |] are easy 
to handle, a singularity like the one that develops in Fig. |^ would soon limit the accessible 
values of p, and thus prevent us from reaching asymptotically high Reynolds numbers. 

At this point, an important difference between the full three-dimensional plane Couette 
problem and the restricted, two-dimensional problem without spanwise degrees of freedom 
shows up. The cumbersome denominator <f>', which does not appear in the restricted case, 
emerges when one eliminates v x (z) from Eq. flT9|) in order to derive the Sturm-Liouville 
equation (p4[) . In this way one arrives at Eq. (|23|) , whereas in the restricted case (k y = 0) 
the condition (|2H| ) of divergence-freeness yields an immediate connection between v x and 
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v' z . Here lies the reason why an asymptotic theory of upper bounds for the full Couette 
problem, paralleling the one developed for the restricted case in Ref. |[22|| , seems to be out of 
reach. Nevertheless, the model problem has led to the identification of the test profiles (|33"D 
as profiles which will provide the best possible bounds on c e in the limit of large Reynolds 
numbers. 

In Fig. [| we visualize a generic representative of the profile family (|33|) . The profile 
functional fll3"D, which becomes a function of S, p, and n when applied to the class (|33"D, now 
reads 



D{cj>} 



n 2 



2 (2n - 1) 5 



By constructing the lower envelope of all graphs of the type shown in Fig. [I], that are 
produced in the (Re, c e )-plane by the individual test profiles, we are able to compute the 
optimal upper bound on c e obtainable from the profiles (|33"D for Reynolds numbers up to 
Re ^ 10 6 . Moreover, from the optimized profiles we can extract the asymptotic scaling of 
the profile parameters. 

In Fig. |5] we report our findings for the optimized upper bound c^(Re) on c e (Re). 
Remarkably, this variational bound separates with a sharp kink from the lower bound 
£e_(Re ) = Re~ x at the energy stability limit, which is located at 

Re ES = 82.65014887(1), c e (Re ES ) = Re^g = 0.012099191758(2). (35) 

(The numbers in brackets denote the uncertainty of the last digit.) Up to Re E s both bounds 
coincide, see Eq. (0). The kink of the upper bound at Re E s is no artifact caused by our 
particular choice of test profiles ([33|): even though some numerical improvement of the 
upper bound is possible in the regime of intermediate Re (see below), one will still obtain 
a bound that increases for Reynolds numbers slightly above Re E s- This behavior is related 



to the singularity of the balance parameter a at Re E s [|20i[22|l , and finds its explanation in 
a change of the role of the optimal background flow. Below Re E s the optimal background 
flow coincides with the laminar flow, that is, with a stationary solution to the Navier-Stokes 
equations, but above Re E s this is no longer the case. Therefore, the deviations v(x, t) from 
the background flow below and above Re E s have a different character. For Re < Re E s these 
deviations are fluctuations around the physically realized laminar flow. Since they will decay 
at least exponentially in time [23|,24], they yield no contribution to e. For Re > Re E s the 



deviations can no longer be regarded as "fluctuations" in the usual sense, and do contribute 
to e. 

The apparent slight non-smoothness of our bound ~c~ e at Re ~ 200, on the other hand, 
indeed is an artifact; it is caused by the restriction of n to values no less than four. We 
have found empirically that the component k x of the wavevector minimizing i? o {0}(k) always 
vanishes. With this knowledge, the system (^) can be simplified: the 20-component complex 
system is converted into a 20-component real one. As a consequence, the parameter value 
n = 3 becomes admissible, which results in a noticeable improvement of the upper bound 
in the intermediate range Re E s < Re < 700, as shown in the inset of Fig. |^. 

For Reynolds numbers around 740 the upper bound exhibits a pronounced minimum, 
and then ascends to the asymptotic value 
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lim c~(Re) = 0.01087(1). (36) 

This value lies slightly above — but within the uncertainty span of — Busse's asymptotic 
result (|[). Thus, our variational bound possesses the same asymptotic -Re°-scaling as the 
bound provided by the Optimum Theory. In Appendix [B| we treat the variational principle 
without taking into account the condition V • V = in Eqs. flUD- This toy problem provides 
an illustrative example for our methods, and allows a direct comparison of variational results 
obtained numerically with the help of the profiles ( |3"3"D to the corresponding results derived 
by Howard PJll| in the framework of the Howard-Busse theory. Our numerical value for the 
asymptotic bound pertaining to this simplified case, lim^g^oo ~c~ e {Re ) ~ 0.07071, differs from 
Howard's analytical value 9/128 ~ 0.07031 by merely 0.6%, even though our variational 
profiles are specifically adapted to the full problem, and not to this simplified case. 

A key for understanding the overall behavior of the variational upper bound lies in the 
expressions i?o{0}(k), considered for fixed profile as functions of the two-dimensional 
wavevector k. We have already mentioned the observation that for each considered test 
profile the absolute minimum of i?o{0}(k) was attained for k x = 0. This finding may seem 
surprising on the first glance, but becomes at least plausible with the help of the following 
argument: if one neglects the condition V ■ V = in Eq. (0), the minimizing k- vector is 
the zero-vector, as shown in Appendix [B]. Taking into account V • V = then enforces 
that the minimizing k becomes non-zero. But since there is no characteristic length-scale in 
streamwise direction, the minimizing k x remains zero. 

There is a further feature that strongly influences the .Re-dependence of the optimized 
profile parameters, and the scaling behavior of the resulting upper bound ~c~ e {Re ) : the single 
minimizing wavenumber k y that characterizes the variational solution at comparatively low 
Re bifurcates to a pair of minimizing wavenumbers, k Vt i and k Vt 2, at the Reynolds number 

Re B « 460. (37) 

The variational principle possesses the intriguing property that above Res both correspond- 
ing minimal values of i?o{0}(k), R\ and R2, are locked to exactly the same value, i.e., two 
eigenvalues of the eigenvalue problem ( fl^[ ) pass through zero simultaneously . This degen- 
eracy reflects the fact that, physically speaking, the sensitivity of the variational problem 
to the behavior of the profiles in the boundary layers separates from the sensitivity to the 
profiles in the interior. A fairly precise mathematical formulation of this statement can be 



given for the Couette problem without spanwise degrees of freedom . 

The consequences of this bifurcation for the upper bound on c £ are dramatic: its approxi- 
mate Re~ 1 ' ' 4 -dependence that prevails for Re < Res changes to the asymptotic Re -behavior, 
while the bound passes through its global minimum 

- £min ^ 0.00970 at #e mi „« 740. (38) 

In order to convey some feeling for the way the variational principle works, we draw Ro{(p}(k) 
as function of both components of k, for the optimized profiles corresponding to the Reynolds 
numbers ReEs (Fig- § cf. Eq. fl35|)) and Re m i n (Fig. [7[ cf. Eq. (^)), respectively. Figure || 
depicts the minimizing values of k y that belong to the variational bound displayed in Fig. ||. 
Whereas the upper /c y -branch scales proportionally to Re in the limit of large Re, the lower 
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branch approaches a constant value. This value is precisely the one that also corresponds 
to the energy stability limit: 

lim ^ = ^ « 0.4960. 

This finding is rather unexpected, since we have the simple laminar profile at the energy 
stability limit, and seemingly entirely different, non-linear profiles with an almost flat interior 
part in the asymptotic regime. However, the restricted model problem has shown that this 
interior part of the asymptotic profiles can be related to the laminar profile, so that the 
asymptotic identity of k Vj i and Ices can be proven analytically by mapping the corresponding 



eigenvalue problems to each other [22 



It is now of considerable mathematical and physical interest to investigate the Re- 
dependence of the optimized profile parameters. We have to keep in mind that the back- 
ground flow U(x) = t/0(^)x is, by construction, not a horizontally-averaged mean flow. 
On the other hand, the optimal profiles invite comparison with physically realized flow 
profiles [19], even though such a comparison remains speculative. Fig. |9| shows the meta- 



morphosis of the optimal profiles with increasing Reynolds number. By a close inspection 
of the data underlying this figure, we can clearly distinguish four different regimes: 

i. ) For Rces < Re < R&i ~ 160, a deformation of the laminar profile takes place. 

Whereas the thickness 5 of the boundary layers remains fixed at 0.5, and the polynomial 
order n of the boundary segments remains at the smallest possible value 3, the profile's 
slope p in the interior decreases from 1 to almost the limiting value 0. In this regime 
the upper bound q reaches its maximum and starts to descend. 

ii. ) The following regime, Re\ < Re < Re2 ~ 670, is characterized by the development of 

boundary layers. The parameter 5 decreases to its minimal value of about 0.14, with n 
remaining fixed at 3, while the slope p increases again. In this regime the bifurcation 
of the minimizing k y occurs, which alters the previous approximate i?e _1 ^-behavior 
of the upper bound. 

iii. ) In the next regime, Re2 < Re < Re% ~ 1845, we observe a dramatic change of the 

shape of the optimized profile's boundary segments. The parameter n — characterizing 
the shape of these segments by means of their polynomial order — increases from the 
smallest possible value n = 3 to n = 34, thus steepening the profiles in the immediate 
vicinity of each boundary, thereby effectively generating new internal boundary layers 
within the original boundary layers. As a consequence of this structural reorganization 
of the boundary layers, the entire boundary segments, the extensions of which are 
quantified by 5, start to reach deeper and deeper into the bulk and finally join again, 
S taking on its maximal value 0.5. In this regime, the upper bound passes through its 
minimum, see Eq. (^), and then changes the sign of its curvature. 

iv. ) For Re above Re%, we find simple scaling laws for the optimal profile parameters. 

These scaling laws give rise to a power-law change of the profile with increasing Re: 
the slope at the boundary z = (or at z = 1) is given by 

n 

0'(O) ~ — ~ a Re, (39) 
2o 
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while the slope at the midpoint z = \ becomes 

=p~(3Re-\ (40) 

with constants a and (3. Parallel to this power-law change of the profile, the upper 
bound on c e ascends smoothly to its asymptotic value. 

If we had not been forced by practical numerical reasons to constrain the shape parameter 
n to values less than or equal to 1000, we would have preserved the asymptotic i?e 1 -scaling 
of n, while the parameter 5 would have remained fixed at 0.5. In our numerics the optimal 
n becomes 1000 at Re ~ 51880, and keeping n fixed at n — 1000 for even higher Re instead 
of allowing it to increase further then forces 5 to decrease, as follows from Eq. (|39"D . The 
ensuing increase of the asymptotic value of the upper bound ~c~ e (compared to the value that 
would be obtainable if one allowed arbitrarily large n) is quite small, as can be seen from the 
fact that c^(51880) = 0.010832(1) differs merely by about 0.3% from the asymptotic value 



V. CONCLUSIONS 

Figure [TTj gives a synopsis of rigorous upper bounds on the dimensionless rate of energy 
dissipation @ in plane Couette flow that have been found so far. The topmost three lines 
indicate results derived with the help of an over-restrictive profile constraint flIT|,|T8"|,|2T| stem- 



ming from functional estimates. This sharpened constraint over-satisfies the actual spectral 
constraint ([14]) and thus leads to bounds exceeding the asymptotic Busse bound (H) by fac- 



tors of 8.8, 7.9, or 6.6, respectively. Only if the spectral constraint is properly implemented 
and evaluated, as done in this work, one can obtain a rigorous bound that practically co- 
incides with the Busse bound in the asymptotic regime, thereby confirming both Busse's 
pioneering work @,[l(| and Kerswell's recent analysis |IJJ . 



However, it must be clearly recognized that much more has been achieved than the 
confirmation of a 25- year old difficult result. Now a variational principle is available that 
can not only produce asymptotic bounds of high quality, but yields rigorous bounds of the 
same quality in the entire range from low to asymptotically high Reynolds numbers. For 
the plane Couette system discussed in this paper, the bound shows a remarkable structure, 
notably a pronounced minimum followed by a i?e-range between 1000 and 1800 in which 
the bound's curvature changes its sign. This occurs at about those Reynolds numbers where 
typical laboratory shear flows start to become turbulent. It is tempting to speculate about 
this coincidence, but such considerations are beyond the scope of the present work. 

Of considerable interest is the mechanism that determines the variational upper bound. 
As depicted in Fig. [8|, the minimizing wavenumber k y bifurcates at Res ~ 460, giving rise 
to one minimizing fc y -branch that asymptotically approaches a constant, and to another one 
that scales linearly with Re. Since both corresponding minima determine the solutions to 



the variational principle by taking on the same value, as described in Sec. [IV], there are two 



characteristic lengths that enter into the solution, one scaling with Re , the other with Re 1 



The analysis of the restricted Couette problem |22j has led to the conclusion that the first, 



.Re-independent minimum reflects the response of the variational principle to the optimized 
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profiles solely in the interior of the fluid layer, whereas the second minimum reflects the 
response to solely the profiles' boundary segments. Thus, the two lengths entering the 
variational solutions can be identified as basically (i) the extension of the interior flat part 
of the optimized profile, and (ii) the effective width of its boundary segments. It is structural 
insight of this kind that will be required for the formulation of a more refined variational 
principle that might ultimately allow one to decide whether there are corrections to classical 
scaling. 

For the moment being, the comparison of our variational bound with the experimental 
data shown in Fig. 10 demands some modesty. The bound clearly is far from being sharp; 
the difference between the variational bound ~c~ e and the corresponding data measured by 
Reichardt flTSJ for the plane Couette flow or by Lathrop, Fineberg, and Swinney [I3| for the 



small-gap Taylor-Couette system still spans an order of magnitude. In short, one important 
step has been made, but there is still a long way to go. 
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APPENDIX A: COMPOUND MATRIX METHOD FOR THE SIXTH-ORDER 

STURM LIOUVILLE PROBLEM 



Keeping in mind that we are not interested in the calculation of the entire spectrum 
furnished by the Sturm-Liouville eigenvalue problem ( p4[) and (|25|), but in the determination 
of that value of R where the smallest eigenvalue passes through zero, we can set A = right 
from the outset. We thus obtain the new sixth-order Sturm-Liouville problem 



^ - 2 



,(5) 







./ 2 



3k 2 + 



, , (3) J^'V i/) 



+ {4fc 2 ^ + iW'}^ 3 > 



3 k +2k z 



0(3) 



k x 

Y 



fi-f R<p {2) > r 



(2) 



k 6 + k A 
k 



6® 

2 k 4 h i k x R 



(2) ) ; 



i — R 



v z = 0, (Al) 



together with the boundary conditions 

V z( Z )\z=0,l = °> V 'z( Z )\ Z =0,l = °> and 



2 = 0,1 



(A2) 
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The eigenvalues of this problem are those values of R where both Eqs. (|Al]) and (|A2|) are 
satisfied. 

In order to integrate Eq. flA~l]), one has to pose suitable initial conditions at one of the 
boundaries, let us say z = 0. We define the six- vectors 

V z (z) = (v z (z) , v' z (z) , vf\z) - i V M( Z ) , v f) {z) , V M( Z ) , v (s) {z) y , 

and consider the three linearly independent fundamental solutions ~V z ,i(z), V^ 2 (z), and 
^ z,z{ z ) to (|A1|) that emerge from the initial values 

V z ,i(0) = (0,0, 0,1,0, Of, V 2 , 2 (0) = (0,0, 0,0,1, Of, and V z , 3 (0) = (0, 0, 0, 0, 0, if , 

(A3) 



respectively. By construction, each of these three solutions satisfies the conditions (|A2|) at 
z = 0, and the most general solution obeying these one-sided b.c. is just a superposition of 



them. That solution which also satisfies the conditions ( A2|) at the other boundary z — 1 is 
then singled out by the requirement 



Det 



/ «*,2(1) v z , 3 (l) 

<i(l) < 2 (1) <s(l) I o. 

U$(l)- 5F«?1(1) «?2(1)-5F«S(1) «?3(l)-5Sr«?l(l) 



But from the numerical point of view the evaluation of this determinant, which in most 
cases necessitates to subtract large numbers of almost equal magnitude from each other, 
is rather inconvenient. An efficient and numerically stable way to bypass this difficulty 
has become known as compound matrix method in the literature, see, e.g., Ref. |27| for an 
introduction to the method, or Ref. ]22j for applications to the restricted Couette problem. 



For applying this technique to the full Couette problem studied in this paper, we define a 
new vector y(z) = (yi(z) , ... , y 2 o(z)) T , the 20 components of which are proportional to the 
3x3 minors of the 6x3 solution matrix that contains V 2) i as its first, V z ,2 as its second, 
and V 2j 3 as its third column: 

,(2) 1 „,W 



Vi = ^ 9 ^sign(cr) v zM1) v' zM2) 



2 £2 ^.°( 3 ) 



(A4) 



2/20 = J2 signO) 



(3) (4) (5) 
V z,o{l) V z,o{2) V z,(i$)- 



The summations involve all permutations a of 1, 2, 3. The powers of k multiplying each sum 
are determined in the following way: the expression defining y 2 o carries the highest total 
number of derivatives, namely 12; this expression is multiplied by k°. Descending in the 
index from 20 to 1, the number of derivatives is successively diminished; each derivative less 
gives a factor of k more. This guarantees two things: firstly, the initial conditions (|29|) are 
independent of k, and secondly, all components have the same order of magnitude. 



From Eq. ( |Al| ) one then obtains a system of first-order equations that already closely 



resembles the system (ESI) which we have studied numerically, with the only difference that 
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in each component i (i = 1, . . . ,20) the term —3kyi appearing in (|28|) is still missing. It 
is obvious that the initial conditions (^) and the boundary condition (|30| ) follow directly 
from ( |A3|) and the definitions 



This term —3 k i/i in the 2-th component of the system is of particular importance. The 



restricted Couette problem has taught us the lesson [22] that the system of first-order equa- 
tions directly provided by the compound matrix method is not suited for numerical analysis 
in the high-i?e-regime, where one has to monitor large values of k. One rather has to rescale 
each component of the system by a common exponential damping factor. The reason for 
this rescaling is related to the fact that the optimized profiles tend to become flat in the 
interior of the fluid layer when Re becomes high, while the extension of the boundary seg- 
ments, where the profiles have to bend in order to meet the boundary conditions 0(0) = 
and 0(1) = 1, tends to zero. This motivates to neglect the profiles' boundary conditions 
altogether, and to consider Eq. ( |Al| ) with a linear profile <fi(z) = cz, so that 0' becomes 
equal to the constant c, and all higher profile derivatives vanish. In the limit c — * one then 
arrives at 

vf> - 3 k\f + 3 k 4 v^ - k 6 v z = 0, (A5) 

with the b.c. (|A2|) . Applying the compound matrix technique as sketched above to this 
boundary value problem, one finds a 20-component system of first-order equations as in 
the case of a general 0, but now with constant coefficients. This system can be solved 
analytically by standard means, but an accurate numerical solution for large fc's will be 
possible only if the solutions do not grow (or decrease) exponentially. Thus, the intended 
numerical approach forces us to rescale the system such that the largest eigenvalue of the 
rescaled system's coefficient matrix is exactly equal to zero. It turns out to be possible to 
determine the Jordan normal form of the 20 x 20-matrix resulting from Eq. ( |A5| ) analytically; 
the largest eigenvalue is non-degenerate and equals 3k. Hence, the transformation 

y x {z)= yi {z)e-* k *, ... , y 20 (z)=y 20 (z)e- 3kz , (A6) 

has the desired effect. Applying the very same transformation also to the system given by 
the compound matrix method for the case of a general candidate profile (and finally omit- 
ting the tilde-signs), changes that system merely be adding —3kyi to the i-th component, 
and thus produces our system (|28|). We emphasize that the identification of the proper ex- 
ponential scaling factor, and the actual scaling transformation (|A6|) , is crucial for obtaining 



equations that remain numerically stable even in the regime of asymptotically high Reynolds 
numbers. 

A striking feature found in the study of the restricted Couette problem, where the system 



corresponding to (|28|) consists of merely six equations p2| , is the possibility to reduce that 
complex 6-component system to a real 6-component system, thereby halving the number 
of equations. This simplification could be obtained from taking lm(yi(z)) = 0, for all z. 
Even though we have strong numerical evidence that the imaginary part of the component 
yi(z) vanishes also in the present case, we did not succeed in deriving a reduced system 



for arbitrary 0. Hence, we make a virtue of necessity and keep the full system (|28D , but 
exploit our insight for controlling the accuracy of the numerical solutions by monitoring the 
magnitude of Im(y 1 (z)). 
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APPENDIX B: NEGLECTING V V = 



In this appendix we study the variational problem for energy dissipation in plane Couette 
flow without requiring that the solutions V(x) to the eigenvalue problem ([14]) be divergence- 
free. Then the solutions to this eigenvalue problem are determined in an enlarged space of 
functions, so that the bound obtainable for this simplified problem can naturally not be as 



good as the one reported in Sec. |TV|. However, the distinct value of this toy model lies in 
two points. Firstly, we can illustrate, condensed in a nutshell, the principles of both the 
background flow method itself and the use of the compound matrix technique. Secondly, this 
model allows an immediate quantitative comparison of our numerically calculated, asymp- 
totic dissipation bound with the corresponding value provided by the Howard-Busse theory. 

To begin with, the Eqs. (|T7D (p0|) for the spectral constraint's eigenvalue problem simplify 
enormously when the condition V • V = is skipped: 

\v x = -2(d 2 z -k 2 )v x + R4>'v z , (Bl) 

XVy = -2 - Vy, (B2) 

Xv z = -2 (d 2 - k 2 ) v z + R <f>'v x . (B3) 

Equation ( P32|) separates from the others, and we immediately obtain a solution that is 
compatible with the b.c. for each A and k, namely v y (z) = 0. Introducing the new 
functions v(z) = v x (z) — v z (z) and w(z) = v x (z) + v z (z), Eqs. flBl]) and ( |B3D transform into 

\v = -2 (d 2 - k 2 ^Jv - R<f)'v, 
Xw = -2 (d 2 z -k 2 )w + Rcp'w; 

the b.c. read v(0) = w(0) = and v(l) = w(l) = 0. 

We now assume that the profile functions <fi(z) satisfy, besides the conditions posed in 
Eq. (|TTD, also the additional monotony condition 

(p'(z) > for < z < 1. (B4) 

Our test profiles (|33J) obviously comply with this plausible requirement. Hence, if the com- 
ponent w(z) does not vanish identically, the positive-definiteness of —d 2 will enforce A > 
for positive R. Consequently, for A = and R > we have w(z) = 0, and are left with the 
second-order Sturm-Liouville boundary value problem 

v" - (k 2 -\R(f)')v = 0, v(0) = v(l) = 0. (B5) 

In exact analogy to the full problem, the main task posed by the spectral constraint is to 
calculate for given and k > the smallest R- value R{<p}{k) where ( [B5|) is fulfilled, and 
then to minimize over all k in order to extract R c {<f)}. But now the second-order differential 
equation (|B5| ) shows that the minimizing k-value k c {<ft} equals zero for each <p{z). In the 
final step, the optimal upper bound on c £ (Re ) follows from inequality ([16]). 

The laminar profile <fr(z) = z yields the "energy stability limit" for our toy problem, 

Re ES = 2tt 2 w 19.74, c £ (Re ES ) = l/2vr 2 w 0.05066. 
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In order to compute R c {<p} for non-laminar profiles we now rewrite the second-order equa- 



tion ( |B5| ) as a system of first-order equations by defining the vector y(z) = (v(z) ,v'{z)) T . 
In this way we get 

2/i = 2/2, y' 2 = -\R<P'yi, (B6) 
where we have set k = 0; initial conditions at z = and b.c. at z = 1 are given by 

y(0) = (0,l) T , 2/i(l) = 0. (B7) 



This reformulation of the Sturm-Liouville problem (|B5|) corresponds precisely to the refor- 
mulation provided by the compound matrix method. If one compares this simple boundary 
value problem ([B6], |B7|) to the corresponding boundary value problem (|28l) -([3~0"l), one gets 
a vivid impression of the complications introduced into the variational principle by the 
condition V ■ V = 0. 

Figure [11] shows the numerically computed variational upper bound on c e (Re) that 
results from inequality (fL2|) when V • V = is neglected, i.e., when the spectral constraint 
is translated into Eqs. ( |B~6| , |B~7| ), and the previous test profiles fl33|) are employed. For 
simplicity, we have fixed the parameter n to the value n — 2. Interestingly, this bound 
exhibits the same shape as the one that had been found analytically when another over- 
restrictive profile constraint stemming from functional estimates had been used instead of 
the actual spectral constraint; that bound is indicated by the long-dashed line in figure [TC . 



The best possible asymptotic upper bound c £ (Re ) for the toy model without V • V = can 
be found for a non-integer value of n, namely m 3.40: 

lim — £ (Re ) « 0.07071. 

Re— >oo 



This value has to be compared to the bound calculated analytically by Howard |q,[TT|j, also 
without accounting for divergence-freeness: 

lim - £ {Re) = w 0.07031. 

Re— >oo iZo 

The work of Kerswell [19| suggests that both the Howard-Busse theory and the background 



flow method, when fully exhausted, give the same asymptotic upper bound. Even though 
we have deliberately used test profiles which are not specifically adapted to the present toy 
model, we have missed Howard's bound by merely 0.57%. 
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FIGURES 

FIG. 1. Graph of an upper bound on c £ produced according to Eq. ( |l6|) by the profile </> shown 
in the inset. 

FIG. 2. Behavior of 4>^(z) /<j)'{z) in the interval < z < 0.015 for 6 = 0.5 and n = 1000. From 
left to right: graphs that result when p is successively diminished by factors of 1/10 from p = 0.1 
top= 10" 6 . 

FIG. 3. Behavior of (jft>(z) /4>'(z) in the vicinity of z = 6 for 5 = 0.002 and n = 4. From top 
to bottom: graphs that result when p is successively diminished by factors of 1/10 from p = 0.1 to 
p = 10" 6 , cf. Fig. |. 

FIG. 4. Test profile cft(z) as given by Eq. (p3D, for a generic parameter set (n = 4, S = 0.3, 
p = 0.1). 

FIG. 5. Bounds on c £ for the plane Couette flow. Points denote the variational upper bound 
c^(Re ) computed numerically with the test profiles (^) ; the solid line on the left is the lower 
bound c^(Re) = Re -1 . The asymptotic value of the upper bound, \miR e _> 00 ~c^(Re) = 0.01087(1), 
lies slightly above, but within the uncertainty span of Busse's asymptotic result (||). — The inset 
depicts the improvement obtained by preassuming that the wavevector minimizing Ro{(f)}(k) has 
a vanishing component in streamwise direction. Solid line: lower bound c £ (Re); dashed line: same 
upper bound as in the main figure; points: improvement of the bound obtained by admitting the 
parameter value n = 3. 

FIG. 6. Three-dimensional surface plot visualizing the k-dependence of Ro{(f)}(k) for the lam- 
inar profile, which is the optimized profile for Reynolds numbers up to the energy stability limit 
Ress- 

FIG. 7. Three-dimensional surface plot visualizing the k-dependence of i?o{</>}(k) for the opti- 
mized profile belonging to the absolute minimum of the upper bound on c E at Re = Re m i n « 740. 
Note the appearance of two distinct minima on the k y -axis. Since the profile results from the 
variational principle, both minima have the same value. 

FIG. 8. Minimizing wavenumber(s) k y corresponding to the upper bound on c £ displayed in 
Fig. i 

FIG. 9. Three-dimensional plot visualizing the metamorphosis of the optimized variational pro- 
files with the increase of the Reynolds number. We have depicted the most important i?e-interval 
on a logarithmic scale, beginning with the energy stability limit ReES an d ending in the scaling 
regime Re > Re^. 
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FIG. 10. Comparison of bounds on c £ (Re) for the plane Couette flow. - Solid slanted 
straight line: lower bound cJ^Re) = Re -1 . - Topmost horizontal solid line: upper bound 
obtained by Doering and Constantin in Refs. [^4|p^ j with the help of an over-restrictive pro- 
file constraint and piecewise linear profiles; c^(Re) m 0.088 for Re > 11.32. - Dashed-dotted 
line: improved bound on c e (Re) derived by Gebhardt et al. in Ref. [18] from the analytical solu- 
tion of the Doering-Constantin principle with the over-restrictive constraint; c^(Re ) ~ 0.079 for 
Re > 16.98. - Long-dashed line: further improvement due to the introduction of the balance 
parameter in Ref. [20]; c^{Re) — > 0.066. This bound is still calculated analytically utilizing the 
over-restrictive profile constraint. - Heavy dots: upper bound for the plane Couette flow ob- 
tained in this work from the variational principle (|l^) with the actual spectral constraint (|i~4]), cf. 
Fig. ||; c^(Re) — > 0.01087(1). The variational profiles are given by Eq. (|33"|), including n = 3 as 
an admissible parameter. — Joining dashed line: asymptotic upper bound (||) derived by Busse 
in Refs. ||[n|; c^(Re) — > 0.010(1). The shaded area denotes the estimated uncertainty of this 
bound. — Triangles: experimental dissipation rates for the plane Couette flow measured by Rei- 
chardt [12|. — Circles: experimental dissipation rates for the Taylor-Couette system with small 
gap as measured by Lathrop, Fineberg and Swinney [ 13 [ . 



FIG. 11. Bounds on c e (Re) for the plane Couette flow obtained by neglecting the condition 
V- V = in the eigenvalue equations (14). Points denote the variational upper bound c^(Re), 
computed numerically from the test profiles ( |33| ) with fixed parameter n = 2; the solid line on the 
left is the lower bound c £ (Re ) = Re^ 1 . 
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